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Abstract 

We report for the first time exact solutions of a completely in- 
tegrable nonlinear lattice system for which the dynamical variables 
satisfy a g-deformed Lie algebra - the Lie-Poisson algebra sug{2). The 
system considered is a g-deformed lattice for which in continuum limit 
the equations of motion become the envelope Maxwell-Bloch (or SIT) 
equations describing the resonant interaction of light with a nonlinear 
dielectric. Thus the A^-soliton solutions we here report are the natural 
g-deformations, necessary for a lattice, of the well-known multi-soliton 
and breather solutions of self-induced transparency (SIT). The method 
we use to find these solutions is a generalization of the Darboux- 
Backlund dressing method. The extension of these results to quantum 
solitons is sketched. 
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1 Introduction 



The Maxwell-Bloch (MB) system of equations have been fundamental to 
much of theoretical quantum optics and nonlinear optics since they were first 
introduced in the late 1960's (some history is in and also in 0). These MB 
systems are of abiding theoretical interest. A feature is that their complete 
integrability is handed down from the 'reduced MB' or (RMB) equations to 
the envelope MB (or self-induced transparency (SIT)) equations, thence, at 
resonance, to the Sine-Gordon equation (cf. e.g. Q). Each member of 
this hierarchy has important physical applications while the physics of SIT 
in particular remains a very active field of current research , even into the 
femto-second pulse regime ||. 

Our recent paper followed up ideas of quantum groups and their rele- 
vance to integrable systems theory and derived a g-deformed lattice version 
of the envelope MB system together with its zero-curvature representation: 
in continuum limit these lattice equations become the resonant envelope MB 
(or SIT) equations. In this paper we now report exact A^-soliton solutions 
of this g-deformed dynamical system. Solitons of the lattice equations were 
promised in and a Riemann-Hilbert method of solution sketched. But for 
the pure A^-soliton solution reported in this paper it is more convenient to 
use a variant of the Darboux-Backlund dressing method which (see below) 
makes an ansatz for the dressing in terms of appropriate N bound states 
eigenvalues. 

Historically P, multi-soliton solutions of the SIT equations were found by 
the method which become Hirota's method 0; Lamb |^ gave the inverse 
scattering solutions; the inverse method for the RMB equations was in 



and obtains the multi-soliton solutions for the SIT equations from it; []T0 
gave a further account of inverse scattering for these SIT equations. These 
several results on inverse scattering confirmed the generality of a method first 
devised to solve the Korteweg-de Vries equation 0. 

Expressed in terms of the complex slowly varying envelopes for the electric 
field and polarisation e, p and the real inversion J\f, the SIT equations can 
be put in the form (e.g. [|10|] ). 



d^e = (p), 

drp + 2i7]p = J\fe, (1) 
drM = -l{e*p + ep*). 
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Propagation is along a coordinate z and ^ = i7x with x = z/c. The time r 
is a retarded time, r = — x); 77 = (cu — uJo)/2Q is the detuning, and fl = 
2TrnoUJofi'^ /h is the couphng constant. The number rio is the density of 2-level 
atoms with the non-degenerate transition frequency uq (rad ■ sec~^); /i is the 
matrix element for dipole allowed transitions at ujq. The star denotes complex 
conjugation and (.) = I^oo^iv){-) d-V is the average over inhomogeneous 
broadening: h{ri) is a 5-function in the sharp-line limit case 0]. 
In m we constructed the completely integrable lattice system whose equa- 
tions of motion for three dynamical variables s„, Hn and (3n at each lattice 
site n can be written 



dtsn = -^(/5„ + /?.+i)(l + 27s„<) + ig2(Mn)(^^ + ,^_^) (2) 

Here g^^" = 1 + 27/3*/3„, q = e"' and 7 > 0, is a real parameter (a coupling 
constant - see below). Reference to ^ shows that in Eqs. (||) we use s„ = 
\/27Xn + i(l'^^"Pn- in [0] the second equation is for dtXn- As can be checked 
(and cf. when the lattice spacing A — for a continuum limit with 

t ^ tA-\ X = nA, (3n = VAS{x), Xn = AS{x), = AS\x) , . 

7 = kA/2, k>0. ^ ' 

one reaches the resonant sharp-line form of the envelope MB (or SIT) equa- 
tions Eqs ([H) via the definitions 

£(e, r) = 2£ix, t), p(e, r) = -2^5(x, t), iV(e, r) = 2S%x, t), (4) 

with = -y//?. Our use of 'lattice Maxwell-Bloch system(LMB) equations' 

for Eqs. (^ stems from this fact. 

A Hamiltonian for this LMB system is 



+ iq 



1 

2 ^ 

^ n=l 

2Hn 



27 X:(/5n+l + Pn) + XniP*n+l + Pn) 
n+lPn - Pn+lPn) } • 



(5) 
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For M < oo it would be natural to impose periodic boundary conditions. 
But we shall look for lattice soliton solutions and here think of M ^ oo with 
suitable boundary conditions still to be specified. The Poisson brackets of 
Eq. (H) are 

{^n?^m} — '>'{'^Hn}6mn, {Hn, Xm} = — i^n^mn, (6) 

and the quantities X*, X„ and form the sUq{2) Lie-Poisson algebra for 
by {•} we mean {x} = {q^ — q~^)/{2^). This algebra has a central element 
X^X* + {Hn}^ = {S}\ The variables X;, X„ enter (D via Xn = [Q]. 
The variables /3* (the 'electric fields', see Eqs. above) satisfy the 

Lie-Poisson g-boson algebra 

(3*^} = tq^^'-d^n, {iVn, (3m] = -tf^n^mn- (7) 



2 The (/-deformed solitons 

In [Jl| we obtained the zero-curvature representation of the system (|^) which 
means that we constructed an over-determined linear system for a matrix- 
function \E'n(C,t) such that 



= L(C|n)^„ (8) 

dt^n = V{C\n)mn, (9) 

where 

2 „-Afn--f/„ 2 

^(Cl^) = E C^'^.H' ^(Cl'^) = — E C'Hn)- (10) 

Here 

V,{n) = 2z7(/5„<_i + I3lsn-i)a\ V^2 = t\<t\ (H) 



while 

A 



Lo(n) = 2z7 ( ° ) - L^^ = \{o^ ± /), (13) 
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The parameter ( & C which appears in Eqs. (|8|)- ([10|) will be thought of as 
the spectral parameter while in continuum limit is a spectral problem in 
L in the usual 2x2 sense (Zakharov-Shabat linear system 0); a^'^'^ are 
the Pauli matrices. The compatibility condition of the two linear systems 
Eqs. (H),(0) under isospectral condition dtC = is 

dtL{C\n) + L{C\n)V{C\n) - V{C\n + l)L{C\n) = (15) 

and this coincides with Eqs. (0), independent of (. However, ( = e^'^^, A G C 
as it was introduced in [0; A is a second 'spectral parameter' and the real axis 
in the A-plane is the circle of the unit radius in the C-plane; A is the usual 
spectral parameter for the equations Eqs. (|l|) derived in continuum limit. 
Notice that time t is suppressed in Eqs. (H),®: an explicit time dependence 
will be indicated only where and when it is needed. Reference to Eqs. (§) and 
(H) may make plain that the function \I/n(C) possesses essential singularities 
of the rank 2 at C = 0, oo. It is also important to notice that the linear 
equations Eqs. (|^) and (^ are invariant under the transformations 

vi/„(C) (-1)"-V^v[;; j^l^ vl/„(C) ^ a^^r.{-Oa^ (16) 

We can now turn to the derivation of exact solutions of the LMB sys- 
tem Eqs. @). For this, as mentioned, we develop a variant of the 



Darboux-Backlund dressing procedure[ll] rather then any inverse scattering 



method [H, |T^. The essence of the dressing procedure is to choose a 'seed' 
solution of the system Eqs. (^, typically some trivial solution, and construct 
from it a new solution associated with additional points C,^, v = 1, . . . ,N 
(say) of the discrete spectrum: thus det 'ifn{Cu,t) = 0, |n|, ^ for the new 
solution \E'„(C, ^)- 

For initial and boundary conditions observe that for SIT and envelope MB 
system Eqs. (P, typical experimental situation is the half-space problem: an 
initial optical pulse enters, supposedly without reflection from x < into the 
resonant medium x > and here breaks up into background radiation and a 
sequence of soliton pulses. The corresponding mathematical problem is the 
Cauchy problem at the point x = 0: e{x,t)\x=o = £oit) together with the 
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asymptotic boundary conditions (in t) that for x > A/" — >■ M-, p — > as 
t — oo. For the so-called 'attenuator' is the ground state iV_ = —1 of 
the inversion density. For the lattice problem we therefore take the half-space 
problem in which /3„(t) and s„(t) are sufficiently decreasing for \t\ — * oo, while 
Hn{t) H such that H corresponds to A^_. In this way we would look for 
a solution in the half-space n > 0, for which it becomes the Cauchy problem 
specified by the conditions 

Pn{t)\n=l = /?l(t), Sn{t)\n=l = Sl(t), Hn{t)\n=l = Hi{t). (17) 

With this as motivation we report in this paper exact A^-soliton solutions 
derived by the dressing procedure based on the seed solution 

f3n = 0, s„ = 0, = H. (18) 

the corresponding solution of the linear system Eqs (||),(|^) is then 



where z{C,) = ^ [C^q ^ — q^^, while the corresponding operator V^^^{(\n,t) 

has = Vi"^ = 0, and V^^ = tW. For the A^-soliton solution of Eqs. (g) 
we construct the new solution ^^-^'{Q of Eqs. (||),(P) through the ansatz 

^fnC) = i^(C)^i°^(C). (20) 

The function F{() is to have poles only at the essential singularities of \l/n(C)- 
As was indicated above these points are and oo. This suggests the ansatz 

AT 

F(C, n, t) = Fo{n, t) + CF+i{n, t) + C'F_,{n, t). (21) 

i=l 

It is convenient to impose the additional conditions on F{C,) that 

^yp* (^1^ ay = F(C), a^F{-Oa^ = (-l)^F(C) (22) 

which are obviously compatible with the transformation Eq. ([16|) We can also 
normalize the matrix F{Q so that for each [n, t) 
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_^(f{n,t) A ( f-\n,t) \ 

^-^-^[ /-(n,t) I, f{n,t))^ ^23) 

where the constant Q is independent of n and t and /(n, t) is a real function. 
We now choose a set of points {Cu}^=i where det ^I/i^^(C) is to vanish. This 
means 



FiCuMCu) = F (^1^ a^$*(C.) = 0, 



(24) 



where 

and constants independent on n and t. The set {Cu,Cu}^=i together 

constitute a necessary and complete set of parameters (spectral data) for an 
A^-soliton solution 0. 

The system of equations Eqs. (^4]) has a unique solution satisfying conditions 
Eqs. (i2D,(i3D if we choose 

Q = gifAf+iEr=i«-, (^ = 6^"-+'"^ (25) 

where •y,^, G R. In so far as C = ^^"'^ = e'^""'"*'^'' and A is the spectral 
parameter for Eqs. (|1|) we are interested in zeros Ci^ defined by the half A- 
plane ImA > which lie inside the circle \C,\ = 1 in the C-plane. The linear 
system Eqs. (^),(|^) is invariant under the gauge transformation Eqs. (^) 
with the potentials written as 

We turn next to a determination of the matrices F±i{n,t). The conditions 
Eqs. (|2^ ) suggest that we should take the matrices F_Ar+2fc diagonal, and the 
matrices F_Ar+2fc-i off-diagonal in agreement with the first of Eqs. (p6|). 
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So will be diagonal from Eq. ( pSj) we set 

FrAi^-7V+2.-l = ? ) , F_-A^-7V+2. = ^q' ) (27) 

in which y^, y^, a;^, are (so far) arbitrary independent complex numbers. 
Then the conditions for the zeros Eqs. ( plf ) mean that we can instead solve 

N 

{l,0) + Y.z,al' = 0,u = l,..,N- (28) 

k=l 

in which the Zk are row- vectors Zk = {xk, yk), and the matrices cr,^ = S^A^S^^ 
in which A^, = diag(C,^, 1/C); the matrices Si, are defined as 



- 1 f $2(0 C$t(C.^ ' ^^^^ 



'^'1(0 -$2(C.. 

and are determined from Eqs. (P^ using the seed solution Eq.([T^). In this 
way the A^-soliton solution of Eqs. (^ is put in the form 

n ^ * 1 Q* 1/1* .on^ 

Pn = -^VnXn, = (30) 



^-2(Af„+//„) ^ ^-2H ^N{n + 1) 
XAr(n) 



where XN{n,t) etc. is determined from Eqs. ( p8D 

For the one-soliton case, = 1, we can choose the single point of the discrete 
spectrum = e'^°~^^°'° (say) and 70 < 0. We then find the formulae 

^ • /2" . , /o ^ expz(0(n,t) - Qq) 

u\ /2" • i,/o xexpi(0(n,t) + ao) 
'-'^'^ = -V^^^"^(^^°) eosh(^(n,t)+7o) ^''^ 

^2{N„+H„) ^ ^2g l - tanh (V;(n,t) - 70) tanh^p 
1 — tanh ('?/'(n, t) + 70) tanh -do 



Here 



[n 



t) = t cosh (270) sin (2ao) — uqq + 0o, (35) 
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ipi^n, t) = t sinh (270) cos (2q;o) — ri'do + '?/'o, (36) 



i?o = 
^•o = 

<^o and ipQ are arbitrary real constants. 

The formulae for various multisoliton solutions for the lattice system Eqs. (|^) 
are too complicated to be presented in detail here. Since for the lattice 
these solutions depend explicitly on the deformation parameter q = e"^, these 
A^-soliton solutions (A^ = 1,2,...) are naturally thought of as q-deformed 
solitons. In continuum limit Eq. (^), g ^ 1 and 7 — and g-deformation 
disappear. 

3 Conclusions and discussion 

As was mentioned above in the case of the real (imaginary) dynamical vari- 
ables and in the sharp line limit the MB system Eqs. ([1|) is equivalent to 
the Sine-Gordon equation. The same procedure is applicable to the LMB 
system which means that in the case of the reduction to the real (imaginary) 
dynamical variables the LMB system is in fact a new version of the lattice 
Sine-Gordon equation. The dressing procedure described in this paper can be 
extended to this case delivering the whole variety of solutions of the (lattice) 
S-G equation (solitons, breathers etc.). In so far as Eqs. (P2D-(PHD form a q- 
deformed soliton we can use this result to gain insight into the quantum case. 
One objective of the investigation of the quantum MB system must be to find 
out, the precise nature of, and to calculate, the 'quantum soliton' solutions. 
In Refs. |l|],[|l4| we introduced and solved exactly through the quantum in- 
verse method (up to the solutions of the Bethe equations) a quantum version 
of the LMB system Eqs. (|). Since this model provides a natural and exactly 
solvable lattice regularization of the continuous limit quantum envelope MB 
(or SIT) system (and recall that the quantum Sine-Gordon can be embedded 



1 sinh (70 - i?7) + sin^ao , ^ 
- In 27 ^ 2 ^ ^7o 

2 sinh (70 + H'^) + sin 

sinh (70 - H-i + iao) 
smh (70 + -H7 + zao) 



(37) 
(38) 
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in this quantum MB) it is very useful for the construction of the evolution op- 
erator and for investigating the quantum dynamics of these continuous mod- 
els which have direct physical meaning. It is known from a number of quan- 
tum models [1^, |16], [l^ that a 'string solution' of the Bethe equations for the 
quantum model corresponds, in the limit of a large number of collective exci- 
tations M, to the soliton solution of the classical counterpart of the exactly 
solvable quantum system. A plausible conjecture which we will justify else- 
where is that the soliton solution Eq. (p2D for the 'electric field' is given by the 
matrix element limM^oo(0|C(Ai)C(A2)...C(AA/)/9i5(Ai)fi(Ai)...S(AM-i)|0), 
where -B(A) is a creation operator for a quasiparticle and C(A) = 5^ (A) is an 
annihilation operator. The rapidities {Xi}fli are roots of the Bethe equations 



M 



2nMA„ sm 7(A.-^g) _ sm7(A,-A,-i) operator (3l is the electric field 

operator which satisfies the g-deformed g-boson algebra analogous to the al- 
gebra Eq. (^. In Ref. ||T7[| it was shown that the creation operator -B(A) plays 
the role of the quantum counterpart of a Blaschke multiplier which builds 
the classical soliton solution |12|. The dressing operator F{Q Eq. (pT]) up to 
certain modifications has the same meaning. This indicates very well how the 
experience obtained in the analysis of the c-number system reported in this 
paper can be used in understanding of the quantum case. In practice this 
experience helps us to trace out the formation of the classical optical soliton 
from a large number of quantum collective excitations, a physical problem of 
considerable interest HTSi . 
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